% initialize globals:
% omega, h, k...


global v0;
v0 = 2;

global v1;
v1 = 2;

global v2;
v2 = 1;

global k;
k = 4;

global omega;
omega = 0.8;

global h;
h = 0.5;

global A;
A = cell(k);

global P;
P = cell(k);

global dimm;
dimm = cell(k);

global D;
D = cell(k);

global iD;
iD = cell(k);

global C;
C = cell(k);

